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Abstract. - We use video microscopy to study a two-dimensional (2D) model fluid of charged 
colloidal particles suspended in water and compute the pressure from the measured particle 
configurations. Direct experimental control over the particle density by means of optical tweez- 
ers allows the precise measurement of pressure as a function of density. We compare our data 
with theoretical predictions for the equation of state, the pair-correlation function and the 
compressibility of a hard-disc fluid and find good agreement, both for the fluid and the solid 
phase. In particular the location of the transition point agrees well with results from Monte 
Carlo simulations. 



Hard-disc (HD) fluids play a prominent role in liquid state theories. This is due to the fact 
that, first, they often serve as reference systems in perturbation theories of two-dimensional 
(2D) liquids (just as hard-sphere fluids do for liquids in three dimensions), and that, secondly, 
at high densities the behavior of every 2D fluid is dominated by excluded volume effects, 
which in turn depends just on the short-ranged hard-core part of the interparticle potential. 
Mainly for theses two reasons, the HD equation of state (EOS) appears also in many theories 
on monolayer adsorption on solid surfaces 0Gj, an aspect illustrated for example in Ref. [2j 
where the HD EOS is used in statistical mechanical theories modeling the binding of peripheral 
globular proteins on lipid membranes. The important role of the HD system explains the 
overwhelming number of theoretical studies on the EOS of a HD fluid, starting as early as 
1959 with the presentation of results from scaled particle theory [3]. Most of all approaches 
to the EOS are based on particular resummations of the virial series and the construction of 
sophisticated Fade approximants [516117) . others use results from integral equation theories [Hj, 
or from theories based on overlap volume functions 9 . Not only the EOS, but also the 
structure of a hard-disc fluid has been explored in detail and is to date well understood 
|51l7ll9l[Tn] . Hard discs are popular also as model system to test advanced density functional 
theories fTTHT^] . used, for instance, to study discs in cavities J3]> or induced freezing and re- 
entrant melting ^3]. Many computer simulation studies of HD systems are available |U15llfi| . 
of which the more recent ones have focused mainly on the melting properties of HD solids |lf>| . 

All these theoretical efforts contrast with the situation on the experimental side. To our 
knowledge, the HD EOS has never been tested with experimental data. The present Letter 
aims at closing this gap. We report on experiments with a 2D model liquid of charged colloidal 
particles suspended in water. Correlation functions computed from real-space images, together 
with the virial equation, are used to calculate the pressure p of the 2D liquid as a function 
of the 2D particle density p, which in our experiment can be conveniently varied by means 
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of optical tweezers. We have thus realized a 2D colloidal model fluid for which the EOS, i.e. 
the p( j o)-diagram, can be determined directly. Comparing the experimental to the theoretical 
EOS of HD's, we find the colloidal liquid to behave as a 2D fluid of HD's over a wide density 
range from the fluid to the solid phase. 

We should remark that monolayers of atoms adsorbed on solid surfaces behave in certain 
cases also as 2D fluids In order to obtain experimental EOS's in these atomic systems 

one has to transform measured isotherms to p(p) diagrams. This has been done by Glandt 
et al. ^7] who compared various theoretical EOS for a 2D Lennard-Jones (LJ) fluid with 
the 2D pressure of argon and krypton, adsorbed on graphitized carbon black ^H]. Still, 
treating adsorbed monolayers as 2D fluids remains an idealization which is rarely justified 
because the 2D system is usually not isotropic due to the natural corrugation of the substrate 
surface EldE]- 

Contrary to atomic systems, where one is simply stuck with the interaction dictated by the 
electronic structure of the atoms, the interactions between colloidal particles can be externally 
controlled and thus adapted to the problem one wants to study. To experimentally realize 
a hard disc system with colloidal particles it is important to assure that the inter-particle 
potential is extremely short ranged and that no attractive parts in the potential exists. We 
have decided not to work with sterically stabilized colloidal particles as this stabilization 
usually leads to a structured pair potential in the distance region were the polymer brushes 
start to overlap [201 • Rather we chose to use charge-stabilized colloids at moderate salt 
concentrations. A screening length kT 1 between 100 nm and 200 nm yielded optimal results. 
Higher salt concentrations resulted in an increased unfavorable particle substrate interaction 
and a pronounced attractive van-der Waals contribution in the interparticle potential. For 
the same reason it proved unpractical to use very thin sample cells, as this is known to also 
induce an attractive part in the pair potential |21j . The colloidal system employed consisted 
of charged sulphate-terminated polystyrene spheres of ao = 3 fim diameter (IDC Cooperation) 
and of charged sulphate-terminated silica particles of 2.4 pm diameter (polydispersity < 4%). 
The suspension was injected into a sample cell made of fused silica plates with 200 pm spacing. 
We performed measurements at high and at low salt concentration, i.e. at kcto = 20.0 and 
kcf o = 6.9, using the 2.4 pm spheres in the high-salt and the 3 pm spheres in the low-salt 
measurement. In order to confine the particles to two dimensions the widened beam of a 
frequency-double Nd:YV04 laser was directed from above into the sample cell, exerting a 
vertical light pressure on the particles and pushing them towards the bottom plate of the 
sample cell. The light pressure and gravity on the one side and the electrostatic particle- 
wall repulsion on the other side create a sharp potential minimum in the vertical direction 
in which the particles are trapped and thus effectively confined to two dimensions. We have 
varied the intensity of the laser beam over a wide range and observed no influence on the 
particle structure. Therefore light induced effects on the pair interaction can be ruled out. In 
contrast to the above-mentioned 2D LJ fluids, realized by adsorbed atoms on solid surfaces, 
we here do not need to make special assumptions concerning the properties of the substrate 
surface. It can be safely ruled out in our experiment that the particles' mobility in 2D 
is in any way hindered by a possible periodicity of the substrate surface. Our system is 
almost perfectly two-dimensional with vertical out-of-plane fluctuations of less than 250 nm. 
An important component of our set-up which is indispensable for measuring p(p) diagrams, 
is our method to vary the particle density p. This was achieved by a scanned optical laser 
tweezers briefly described in the following. The beam of a laser was reflected from a computer- 
controlled system of mirrors and focused into the 2D system plane. By means of these mirrors 
the tweezers was repeatedly (200Hz) scanning a line having the form of a rectangular box, 
which results in an optical line trap for the particles along the contour of the box. Trapped 
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particles arrange like a pearl-necklace along this contour and, due to the repulsive inter-particle 
potential, prevent other particles from escaping the box, thus defining a system (all particles 
in the box) with a certain number of particles (1000 to 3000). The box size can be changed via 
the control unit of the mirror system which allows precise adjustment of the particle density 
in the system. For each system at density p, digital video microscopy measurements were 
made using a high-aperture long-distance microscope objective (Zeiss Achroplan 63x) with 
dark field illumination. The particle positions were determined with an accuracy of better 
than 50 nm. From statistical averages of the particle positions, we computed pair correlation 
functions g(r). Further details about the experimental set-up and the data analysis can be 
found in our previous papers |22U23| . 

Having the pair-correlation function at density p, the pressure can be computed from the 
virial equation which for a isotropic 2D system reads 

(ipa 2 = pa 2 - [ drr 2 (3u'(r)g(r) . (1) 

2 Jo 

where (3 = 1/kT. u'(r) is the derivative of the interparticle potential u(r) with respect to 
the distance r; this potential is not known a priori, but may be derived from the structural 
data, e.g. from inverting g(r). This was done by means of the Ornstein Zernike equation and 
appropriate closure relations (Percus-Yevick and HNC), a method described and tested in |22|. 
see also Performing then this integral for all measured g(r)'s, we obtain p(p), the desired 
EOS. For the nag = 20.0 measurement, we obtained from the inversion almost perfect Yukawa 
potentials (u(r) ~ er Kr jr) which were - within the error bars specified below - identical, 
regardless at what density we analyzed the pair-correlation function, while for the low salt 
measurement we obtained, for all g(r)'s, Yukawa potentials with almost identical prefactors 
and screening constants k, but with a density-dependent truncation at large distances and 
high densities. This density-dependent truncation has been interpreted in terms of many-body 
interactions |23II22| . The fact that in both measurements the k resulting from the inversion 
procedure do not show any dependence on p, clearly indicates that the salt-ions dominate the 
screening behavior and that the contribution of counter-ions to the screening can be neglected. 
Due to the double- layers around the particles, the effective hard-core diameter a is larger than 
the actual particle diameter ctq. Similar to |2fi| . we determined a by first evaluating the second 
virial coefficient 

/■OO 

B 2 = -n r(e _ " u M - l)dr (2) 
Jo 

using the interparticle potential u(r) obtained from the measured <?(r)'s, and by then identi- 
fying £?2 with the second virial coefficient of a pure HD system, B^ 10 ^ — na 2 /2. We obtain 
a = 1.084<7o for the kctq = 20.0 measurement, and a = 2.16<7o for the ko-q = 6.9 measurement. 

Fig. Q shows the pressure as a function of pa 2 , for both measurements. The precision 
of our procedure is limited mainly by errors made in the determination of u(r); they lead to 
small variations in the computed pressure, but also in the effective hard-core radii. Therefore, 
the error bars appear tilted in Fig. JTJ. The measured EOS is compared with Monte-Carlo 
(MC) data from and with an expression proposed by Baus and Colot jj] for the EOS of 
a fluid HD system, 

^ = 1+ (Hy" (3) 

with n = npa 2 /A and c x = 0, c 2 = 0.128, c 3 = 0.0018, c 4 = -0.0507, c 5 = -0.0533, c 6 = 
—0.0410. The agreement between our high-salt data and the HD EOS is almost perfect, both 
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Fig. 1 - Equation of state of a hard-disc system in the fluid and the solid density range. Theoretical 
predictions for the fluid branch (Baus/Colot |7j) and the solid branch (Velasco/Mederos 1121) are 
compared with experimental data, measured in high-salt (filled circles) and low-salt (filled squares) 
colloidal suspensions. Error-bars (solid short lines attached to the filled symbols) are inclined for 
reasons explained in the text. Monte Carlo data for a Yukawa fluid (dashed line) are provided to 
interpret the kctq = 6.9 measurement. 



on the fluid and the solid side of the phase transition, but also with regard to the location 
of the transition point itself. A much simpler EOS (Sj (c2 = 0.125 and c n = for n = 1, 
n > 2 in eq. Q) agrees equally well and also the density- functional theory by Velasco and 
Mederos ^2] (data not shown here). Surprisingly, even for low salt concentration (kcto = 6.9, 
c/co = 2.16 !), the colloidal system can still be successfully mapped to a HD fluid, at least for 
low densities. For per 2 > 0.7, there are marked deviations from the HD EOS. We performed 
MC simulations to compute pair-correlation functions of a quasi-2D fluid system, using the 
Yukawa part of u{r) of the kcto = 6.9 measurement. The pressure curve, computed with eq. 
from the MC-generated pair-correlation functions, is plotted as dashed line in Fig. QJ. It is 
evident that this line differs, as expected, from the HD EOS, but also from the experimental 
values of the kctq = 6.9 measurement, a difference which hence must be due to deviations of 
the experimental pair potential from the assumed Yukawa form. As pointed out in |22 , 23 
the experimental pair potential shows a density-dependent truncation in u(r), which has been 
ignored in the MC simulation. 

Baus and Colot [7] have suggested a semiempirical expression for the direct correlation 
function of a HD fluid which can be related to g(r) via the Ornstein-Zernicke equation. Fig. 
shows the resulting pair-correlation functions (solid lines) for various densities and compares 
them with our measured g{r). We emphasize that after matching the second virial coefficient 
no free parameter was used. While in the case of high salt concentration excellent agreement 
is found, differences are observed for the nag = 6.9 measurement, mainly for small particle 
separations where the first peak is less pronounced compared to that of the HD fluid, due to 
the softer pair potential. For larger distances, packing effects dominate, and the agreement is 
again good. 

From the recorded particle positions we can also compute the particle number fluctuation 



Brunner et 



5 




Fig. 2 - Comparison between theoretical [J] (solid line) and experimental (filled circles) pair- 
correlation functions for a hard-disc fluid at different densities, experimentally realized with low-salt 
(kcfo — 6.9) and high-salt (kctq = 20.0) colloidal suspensions. 



(AN 2 )/(N) with AN = N — (N), (N) being the mean number of particles. This quantity 
which is related to the isothermal compressibility, xrkTp = (AN 2 )/(N), is plotted as open 
symbols in Fig. Alternatively, we can calculate \t from the density derivative of the 
pressure, \rkTp = {fidp/dp) -1 , which we did in Fig. J3J| (filled symbols) computing the dif- 
ferences between neighboring experimental values in Fig. (JTJ . The compressibilities, calculated 
in both ways, are compared with the derivative of Henderson's EOS 6 

[ 1 + 7778 V + 8 r 1 

pkTXT = [ji^W +T1 4ti^r\ ' (4) 

plotted in Fig. © as solid line. While the density derivatives of the experimental pressure 
values agree nicely with the theoretical prediction of eq. (AN 2 )/(N) shows deviations, 
especially for the low-salt measurement. One reason for this discrepancy might be a finite 
size effect. To estimate this effect, we used the pair-correlation functions, suggested by Baus 
and Colot for the infinite system, to calculate the particle number fluctuation in a finite 
subvolume (V — tttq, tq — lOoo) of a hard-disc fluid composed of a fixed number of particles 
(1000 particles), in a way described in detail by Roman et al. |2Z|. The result is plotted in 
Fig. @ and the small correction illustrates that, at least for the high salt measurement, the 
observed differences can be explained with a finite size effect. For the low salt measurement, 
the remaining differences are probably due to insufficient sampling. 

To obtain more information on the character of the different phases in Fig. JJ), we also 
examined the orientational correlation function ge{r) |24| for the high-salt measurement. Up 
to a density of pa 2 = 0.85 the system exhibits a pure liquid phase as confirmed by the isotropic 
pair-correlation function g(r) and the exponentially decaying ge(r). On the other hand, at 
densities above pa 2 = 0.95 the system is in stable crystalline state, with a slow algebraically 
decaying pair correlation function and a constant orientational correlation function. However, 
in between, i.e. for 0.85 < pa 2 < 0.95, there is a transition region, in which ge(r) decays 
slower than exponentially and a local hexagonal order can be observed. At present it cannot 
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Fig. 3 - Compressibility of a hard-disc fluid, as obtained from density derivatives of the p(p) curves 
of Fig. Q and the particle number fluctuations {AN 2 )/{N}. Symbols for the experimental values, 
lines for theoretical predictions for an infinite system, eq. Q, and a finite system |27) . 



be decided whether this is the hexatic phase or a very polycrystalline state. In contrast to 
the measurement of Marcus and Rice (20j no pronounced liquid-hexatic and hexatic-solid 
coexistence region has been observed, which is most likely due to the fact that their pair- 
potentials had weak attractive parts, while ours are purely repulsive. 

We close with a few comments and summarizing remarks, (i) Similar hard-disc-like fluids 
have been experimentally realized also by Murray et al. [2H| and by Marcus et al. [201 using 
colloidal suspensions in confined geometries. While these authors concentrated on the melting 
properties of 2D systems, it has been the focus of the present work to study the EOS of the 2D 
system which naturally requires a convenient experimental tool to vary the density, realized 
here with optical line tweezers, (ii) To experimentally realize a HD system, it is decisive 
to choose the right salt concentration, as salt tunes the inter-particle potential but also the 
interaction between particle and wall. Too much salt leads to an attraction in the interparticle 
potential and makes the wall-particle interaction more short-ranged (so that the 2D system 
plane is shifted too close to the wall), while too little salt results in a very soft repulsive 
potential. In both cases, the colloidal system can no longer be interpreted as a hard disc 
system. The extent of the deviation between the experimental and the HD system caused by 
a very soft pair potential (koo = 6.9) has been shown above, (iii) 2D colloidal suspensions can 
indeed be considered as a hard-disc system. This applies both to the structure as well as to 
the thermodynamics. The effective HD radius of the colloidal particles is correctly defined by 
simply demanding the second virial coefficient to be equal in both descriptions, (iv) On the 
atomic level, the EOS of 2D fluids can be derived from adsorption isotherms As opposed to 
these experiments, we here have the full structural information, controllable substrate-particle 
interactions and tunable inter-particle potentials. In adsorbed monolayer films |18| the phase 
behavior can be deduced only via averaged quantities and there is a certain ambiguity about 
the emergence of multi-layer adsorption at higher densities. This highlights once again the 
power of colloids as a model system in statistical mechanics. 



Brunner et al.: The equation of state of a hard-disc fluid 



7 



REFERENCES 
[1] W.A. Steele, J. Chem. Phys. 65, 5256 (1976). 

[2] W.A. Steele, The Interaction of Gases with Solid Surfaces (Pergamon Press, Oxford, 1974). 
[3] R.C. Chatelier and A. P. Minton, Biophys. J. 71, 2367 (1996); T. Heimburg and D. Marsh, 

Biophys. J. 68, 536 (1995); T. Heimburg, B. Angerstein, and D. Marsh, Biophys. J. 76, 2575 

(1999). 

[4] H. Reiss, H.L. Frisch, and J.L. Lebowitz, J. Chem. Phys. 31, 369 (1959). 

[5] F.H. Ree and W. G. Hoover, J. Chem. Phys. 40, 939 (1964); F.H. Ree and W. G. Hoover, J. 

Chem. Phys. 46, 4181 (1967); K.W. Kratky, Physica 85A, 607 (1976); K.W. Kratky, J. Chem. 

Phys. 69, 2251 (1978); D. Henderson, Mol. Phys. 34, 301 (1977); M.J. Maeso, J.R. Solana, J. 

Amoros, and E. Villar, Mater. Chem. Phys. 30, 39 (1991). 
[6] D. Henderson, Mol. Phys. 30, 971 (1975). 

[7] M. Baus and J.L. Colot, J. Phys. C 19, L643 (1986); M. Baus and J.L. Colot, Phys. Rev. A 36, 

3912 (1987), contains a comprehensive discussion of all theories in 0. 
[8] F. Lado, J. Chem. Phys. 49, 3092 (1968); Y. Uehara, T. Ree, and F.H. Ree, J. Chem. Phys. 70, 

1876 (1979); M.S. Ripoll and C.F. Tejero, Mol. Phys. 85, 423 (1995); J. Sedlbauer, S. Labik, A. 

Malijevsky, and W.R. Smith, Mol. Phys. 85, 61 (1995); L. Verlet and D. Levesque, Mol. Phys. 

46, 969 (1982). 

[9] D.J. Gonzalez, L.E. Gonzalez, and M. Silbert, Mol. Phys. 74, 613 (1991). 
[10] N. Studart, H.V. DaSilveira, and U. DeFreitas, Mod. Phys. Lett. B 10, 393 (1996). 
[11] M.Takamiya and K. Nakanishi, Mol. Phys. 70, 767 (1990); C.F. Tejero and J.A.Cuesta, Phys. 

Rev. E 47, 490 (1993); P. Tarazona, Physica A 306, 243 (2002). 
[12] E. Velasco and L. Mederos, Phys. Rev. B 56, 2432 (1997). 
[13] S.C. Kim, Z.T. Nemeth, M. Heni, and H. Lowen, Mol. Phys. 99, 1875 (2001). 
[14] C. Bechinger, M. Brunner, and P. Leiderer, Phys. Rev. Lett. 86, 930 (2001); L.L. Rasmussen 

and D.W. Oxtoby, J. Phys.: Condens. Matter 14, 12021 (2002). 
[15] B.J. Alder and T.E. Wainwright, Phys. Rev. 127, 359 (1962); N.E. Thompson and B.C. Freasier, 

Mol. Phys. 41, 127 (1980); J.J. Erpenbeck and M. Luban, Phys. Rev. A 32, 2920 (1985). 
[16] J.F. Fernandez, J.J. Alonso, and J. Stankiewicz, Phys. Rev. Lett. 75, 3477 (1996); H. Weber, D. 

Marx, and K. Binder, Phys. Rev. B 51, 14636 (1995); K.Binder, S. Sengupta, and P. Nielaba, 

J.Phys.: Condens. Matter 14, 2323 (2002); S. Sengupta, P. Nielaba, and K. Binder, Phys. Rev. 

E 61, 6294 (2000); A. Jaster, Europhys. Lett. 42, 277 (1998); A. Jaster, Phys. Rev. E 59, 2594 

(1999); M.A. Bates and D. Frenkel, Phys. Rev. E 61, 5223 (2000). 
[17] E.D. Glandt and A.L. Myers, J. Chem. Phys. 70, 4243 (1979). 

[18] A. Thorny and X. Duval, J. Phys. Chem. 67, 1101 (1970); F.A. Putnam and T. Fort, J. Phys. 

Chem. 79, 459 (1975). 
[19] PA. Monson, W. Steele, and D. Henderson, J. Chem. Phys. 74, 6431 (1981). 
[20] A.H. Marcus and S.A. Rice, Phys. Rev. Lett. 77, 2577 (1996); A.H. Marcus and S.A. Rice, 

Phys. Rev. E 55, 637 (1997). 
[21] S.H. Behrens, D.G. Grier, Phys. Rev. E 64, 050401 (2001). 

[22] R. Klein, H.H. von Griinberg, M. Brunner, C. Bechinger, V. Lobaskin, J. Phys.: Condens. 

Matter, 14, 7631, (2002). 
[23] M. Brunner, C. Bechinger, W. Strepp, V. Lobaskin, and H.H. von Griinberg, Europhys. Lett., 

58, 926, (2002). 

[24] D.R. Nelson, and B.I. Halperin, Phys. Rev. B, 19, 2457 (1979). 

[25] M. Quesada-Perez, A. Moncho-Jorda, F. Martinez-Lopez and R. Hidalgo- Alvarez, J. Chem. 

Phys. 115, 10897 (2001). 
[26] J.A. Barker and D. Henderson, J. Chem. Phys. 47, 4714 (1967). 

[27] F.L. Roman, J.A. White, A. Gonzalez, and S. Vealsco, J. Chem. Phys. 110, 9821 (1999); we 
used eq.(10) and (32) in F.L. Roman, J.A. White, and S. Velasco, J. Chem. Phys. 107, 4635 
(1997). 

[28] C.A. Murray and D.H. Van Winkle, Phys. Rev. Lett. 58, 1200 (1987). 



